Search for: All records

Abstract The nonlinear Schrödinger (NLS) equation in one dimension is considered in the presence of an intensity-dependent dispersion term. We study bright solitary waves with smooth profiles that extend from the limit where the dependence of the dispersion coefficient on the wave intensity is negligible to the limit where the solitary wave becomes singular due to vanishing dispersion coefficient. We analyse and numerically explore the stability for such smooth solitary waves, showing with the help of numerical approximations that the family of solitary waves becomes unstable in an intermediate region between the two limits, while being stable in both limits. This bistability, which has also been observed in other NLS equations with generalized nonlinearity, brings about interesting dynamical transitions from one stable branch to another stable branch, which are explored in direct numerical simulations of the NLS equation with the intensity-dependent dispersion term. 

IntroductionThe moment quantities associated with the nonlinear Schrödinger equation offer important insights into the evolution dynamics of such dispersive wave partial differential equation (PDE) models. The effective dynamics of the moment quantities are amenable to both analytical and numerical treatments. MethodsIn this paper, we present a data-driven approach associated with the “Sparse Identification of Nonlinear Dynamics” (SINDy) to capture the evolution behaviors of such moment quantities numerically. Results and DiscussionOur method is applied first to some well-known closed systems of ordinary differential equations (ODEs) which describe the evolution dynamics of relevant moment quantities. Our examples are, progressively, of increasing complexity and our findings explore different choices within the SINDy library. We also consider the potential discovery of coordinate transformations that lead to moment system closure. Finally, we extend considerations to settings where a closed analytical form of the moment dynamics is not available. 

Abstract In the present work we revisit the problem of the quantum droplet in atomic Bose–Einstein condensates with an eye towards describing its ground state in the large density, so-called Thomas–Fermi (TF) limit. We consider the problem as being separable into 3 distinct regions: an inner one, where the TF approximation is valid, a sharp transition region where the density abruptly drops towards the (vanishing) background value and an outer region which asymptotes to the background value. We analyze the spatial extent of each of these regions, and develop a systematic effective description of the rapid intermediate transition region. Accordingly, we derive a uniformly valid description of the ground state that is found to accurately match our numerical computations. As an additional application of our considerations, we show that this formulation allows for an analytical approximation of excited states such as the (trapped) dark soliton in the large density limit. 

Abstract The Riemann problem for the discrete conservation law is classified using Whitham modulation theory, a quasi‐continuum approximation, and numerical simulations. A surprisingly elaborate set of solutions to this simple discrete regularization of the inviscid Burgers' equation is obtained. In addition to discrete analogs of well‐known dispersive hydrodynamic solutions—rarefaction waves (RWs) and dispersive shock waves (DSWs)—additional unsteady solution families and finite‐time blowup are observed. Two solution types exhibit no known conservative continuum correlates: (i) a counterpropagating DSW and RW solution separated by a symmetric, stationary shock and (ii) an unsteady shock emitting two counterpropagating periodic wavetrains with the same frequency connected to a partial DSW or an RW. Another class of solutions called traveling DSWs, (iii), consists of a partial DSW connected to a traveling wave comprised of a periodic wavetrain with a rapid transition to a constant. Portions of solutions (ii) and (iii) are interpreted as shock solutions of the Whitham modulation equations. 

Abstract In the present work we revisit the problem of the generalised Korteweg–de Vries equation parametrically, as a function of the relevant nonlinearity exponent, to examine the emergence of blow-up solutions, as traveling waveforms lose their stability past a critical point of the relevant parameterp, here atp = 5. We provide anormal formof the associated collapse dynamics, and illustrate how this captures the collapsing branch bifurcating from the unstable traveling branch. We also systematically characterise the linearisation spectrum of not only the traveling states, but importantly of the emergent collapsing waveforms in the so-called co-exploding frame where these waveforms are identified as stationary states. This spectrum, in addition to two positive real eigenvalues which are shown to be associated with the symmetries of translation and scaling invariance of the original (non-exploding) frame features complex patterns of negative eigenvalues that we also fully characterise. We show that the phenomenology of the latter is significantly affected by the boundary conditions and is far more complicated than in the corresponding symmetric Laplacian case of the nonlinear Schrödinger problem that has recently been explored. In addition, we explore the dynamics of the unstable solitary waves forp > 5 in the co-exploding frame. 

Deflation is an efficient numerical technique for identifying new branches of steady state solutions to nonlinear partial differential equations. Here, we demonstrate how to extend deflation to discover new periodic orbits in nonlinear dynamical lattices. We employ our extension to identify discrete breathers, which are generic exponentially localized, time-periodic solutions of such lattices. We compare different approaches to using deflation for periodic orbits, including ones based on Fourier decomposition of the solution, as well as ones based on the solution’s energy density profile. We demonstrate the ability of the method to obtain a wide variety of multibreather solutions without prior knowledge about their spatial profile.